Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN I^TeX style file v2.2) 



Probing the bispectrum at high redshifts using 21 cm HI 
observations 



o 
o 

(N 

in 

(N 
> 
00 



in 
o 



o 



X 



SK. Saiyad Ali^*, Somnath Bharadwaj^f and Sanjay K. Pandey^J 

^ Department of Physics and Meteorology & Centre for Theoretical Studies , IIT Kharagpur, Pin: 721 302 , India. 

^ Deptt. of Mathematics ,L.B.S. College, Gonda 271001 , India & Visiting Associate, Inter University Centre for Astronomy and Astrophysics, 
Pune , India. 



5 February 2008 



ABSTRACT 

Observations of fluctuations in the redshifted 21 cm radiation from neutral hydrogen 
(HI) are perceived to be an important future probe of the universe at high redshifts. 
Under the assumption that at redshifts z < 6 (Post-Reionization Era) the HI traces 
the underlying dark matter with a possible bias, we investigate the possibility of using 
observations of redshifted 21 cm radiation to detect the bispectrum arising from non- 
linear gravitational clustering and from non-linear bias. We find that the expected 
signal is ^ O.lmJy at 325 MHz [z = 3.4) for the small baselines at the GMRT, the 
strength being a few times larger at higher frequencies (610 MHz, z = 1.3). Further, the 
magnitude of the signal from the bispectrum is predicted to be comparable to that from 
the power spectrum, allowing a detection of both in roughly the same integration time. 
The HI signal is found to be uncorrelated beyond frequency separations of ~ 1.3MHz 
whereas the continuum sources of continuum are expected to be correlated across 
much larger frequencies. This signature can in principle be used to distinguish the HI 
signal from the contamination. We also consider the possibility of using observations 
of the bispectrum to determine the linear and quadratic bias parameters of the HI at 
high redshifts, this having possible implications for theories of galaxy formation. 
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1 INTRODUCTION 

Quantifying non-Gaussian features in the large-scale mat- 
ter distribution is an important issue in cosmology. It is 
widely accepted that the large-scale structures originated 
from initially small density fluctuations which were a Gaus- 
sian random field. Non-Gaussian features arise as the fiuctu- 
atio ns grow throu gh the process of gravitational instability 
Ceg. lPeebleJlQSOri . and quantifying these at various stages of 
the growth is expected to yield a significant amount of infor- 
mation. The three-point correlation function and its Fourier 
transform, the bispectrum, are the lowest order statistics 
sensitive to non-Gaussian features. 

The determination of the three-point corr elation 
function was pioneered bv lPeebles and GrothI lll975ll using 
the galaxies in the Lick and Zwicky angular catalogues. 
Later studies of the th ree point corre la tion function usin g 
galaxy surveys include iPeeblej Jl98 j) . | Bean et al 
lEfstathiou fc Jcdrzcicwskil Jl984i) . iHal^^uttoi^t^l] 
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il989[l . IJing. Mo fc BorneJ Jl99ll). iGaztanaea fc Frieman 
1,199^, iFrieman fc Gaztanagal ^^^^ Ijing fc Borner 
(1998). The recently completed Two -degree Field Galaxies 
Redshift Survey (2dfGRS) and the currently ongoing SDSS 
have been used to measure the galaxy three-point correlation 
function at an unprece dented level of a ccuracy (2d FGRS, 
ping & Borner 2004; Kavo et all l2004l : ICaztanag a et alJ 
2005; SDSS ) allowing various issues like the luminosity and 
colour dependence to be addressed. There have recently 
been a f ew measurements of the bispectrum for the I RAS 



galaxies jFeldman et alJ 



200lt ISco ccimarro et al.ll200lll and 



for the 2dFGRS ( 200 aVerde et a l. 2002) which have (as 
discussed later) been used to estimate the galaxy bias. 

All the measurements of the three-point correlation 
function or the bispectrum mentioned here are restricted 
to low redshifts {z < 1). This limitation arises because they 
are based on galaxy surveys which do not extend to very 
high redshifts. For example the 2dFGRS which is the largest 
completed redshift survey extends to a maximum redshift 
of around z = 0.3. Further, a galaxy or quasar survey ex- 
tending out to large redshifts [z > 1) would select only the 
brightest objects which would be sparse in number. These 
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objects would be distributed in the densest regions which 
are known to have a highly biased distribution and their 
distribution would not be representative of the underlying 
matter distribution. In this paper we discuss how observa- 
tions of the redshifted 21 cm emission from neutral hydrogen 
(HI) can be used to measure the bi-spectrum at high red- 
shifts (ie. z ~ 1 and higher). 

The HI density in the redshift range 1.7 — 5.5 is 
known from observati ons of Ly-a absorption lines along lines 



19961: "Storrie-Lombardi, & Wolf^ 


I2OOOI: iPeroux et alJl200ll 


Prochaska. Hcrbcrt-Fort & Wolfe 




[20051) . These observa- 



tions indicate that Q,g{z), the comoving density of neu- 
tral gas expressed as a fraction of the current criti- 



cal density, has a value 



10" 



Recent observations 



JProchaska. Herbert-Fort fc Wolf3 l2005l) show that ^g{z} 
drops by about 25% at 2 = 2.2 — 2.5. Current observations 
of Q,g [z) are consistent with no evolution at higher redshifts 
(3 < z < 5.5). The bulk (~ 70%) of the HI resides in high 
column-density regions (> 2. x 10^" cm"'^) which are possi- 
bly the progenitors of present day galaxies. The redshifted 
21 cm (1420 MHz) emission from this HI is present as a 
background radiation in radio observations at all frequen- 
cies below 1420 MHz. The fluctuations in this background 
radiation arise from fluctuations in the HI density and from 
peculiar velocities at the red shift where the radiation orig- 
inated y3harad wai. Nath and Sethi 2Q0L hereafter BNS). 
The possibility of detecting this holds the potential of being 
an important tool to probe the large-scale structures in the 
redshift range z = 1 — 6. This is particularly significant in 
the context of the GMRT llSwarup et al.lll99J) which oper- 
ates in several bands in the frequency range of interest. It 
is also interesting to note that the fluctuations in the back- 
ground HI radiation is expected to exceed the fluctuations 
in the CMBR by a factor of 10 — 100 over the frequency 
range of interest (BNS). 

Complex visibilities at different baselines and frequen- 
cies are the primary quantity measured in radio interfer- 
ometric observations. Correlations between pairs of visi- 
bilities directly probe the three dimensional power spec- 
trum of HI fluctuation at the redshift where the radia- 
tion originated. This h as been s tudied in a series of papers 
teharad^^^^^etin[i2001; Bh aradwai and PandevI I2OO3I 
iBharadwai fc Srika nt' 2004') which develop the relation be- 
tween visibility correlations and the HI power spectrum, and 
estimate the expected signal in the redshift range z < 6 
both analytically and using N-body simulati ons. This for- 
malism has later been extended (Bharadwai a nd Ali Il2005l) 
to quantify the expected HI visibility correlation signal 
over a large redshift range extending from the dark ages 
[z ~ 100) to the present epoch. The possibility of using vis- 
ibility correlations to probe the 3D HI power spectrum has 
also been considered by Morales and Hewitt ( 200-1) while 
IZaldarriaea. Furlanetto and HernauistH (j200^ propose the 
use of the angular power spectrum, both in the context of 
det ecting the Epoch of Reionizatio n (EOR) signal. 

IBharadwai and PandevI (|200^ ( hereafter BP05) have, 
in a recent paper, investigated the possibility of using corre- 
lations between visibilities observed at three different base- 
lines and frequencies to measure the bispectrum of HI fluc- 
tuations at the epoch of reionization. The HI distribution 
during the epoch of reionization will be determined more 



by the size, distribution and topology of the ionized regions 
and it is anticipated that this will not trace the underlying 
dark matter distribution. It is thus expected that the bispec- 
trum of the HI during reionization will tell us more about 
the ionization process and less about the non-Gaussianity 
arising from the non-linear growth of density perturbation. 
In this paper we investigate how three visibility correlations 
can be used to study the HI bispectrum at redshifts z < 6, 
where it is reasonable to assumed that the large scales HI 
distribution traces the dark matter with a possible bias. 

We next briefly discuss what we can hope to learn from 
observations of non-Gaussianities in the HI. In perturba- 
tion theory the lowest order at which there is a non-zero 
three point correlation function is the second order 
119841: lBharadwailll994^ . These predictions are expected to 
be valid on large scales which are weakly non-linear, and 
it is in principle possible to compare these with observa- 
tions and test the currently accepted scenario for struc- 
ture formation. The bispectrum has been perceived to be 
a more effective statistics for quantifying non-Gaussianity 
in the weakly non-linear regime and its theory has b een de- 
velope d in .Frv (19941 



. ^ Hiv on ct al. ( 199 5') . Matarrcs ^^^^ 

199aVerde efal] 1998 1998. IScoccimarro et al 



^98 )■ [Scocdnwro ^Cpuchinan fc FriemanI 1^9^ and 



IScocdniarr^^^dT ~ (l20Mll . These investigations show that 
it possible to use observations of the bispectrum (or 
equivalently the three point correlation function) to de- 
termine the linear and quadratic bias parameters. This 
has been actually carried out using various galaxy sur- 
veys like the APM galaxies and the IR AS galaxies 
iFeldman et al. 200 ll: IScoccimarro et al.ll200ll) . In a recent 
analysis 200 z Verde et alj 2002 2002 have analysed the bis- 
pectrum of the 2dFGRS to conclude hi = 1.04 ± 0.11 
for the linear bias parameter and 62 = —0.054 ± 0.08 for 
the quadratic bias parameter, indicating that the 2dFGRS 
galaxies are an unbiased tracer of the underlying dark mat- 
ter. Further, the linear bias parameter was used in combi- 
nation with the redshift distorti on parameter (3 ~ n ^n/fei 
measured from the same survey JPeacock et alJl200lD to de- 
termine the density of the universe fimo = 0.27 ± 0.06. A 
new analysis of t he three point correlati on function of the 
2dFGRS galaxies (ICaztanaea et aLll2005^ confirms that the 
large-scale structures formed through gravitational instabil- 
ity starting from Gaussian initial conditions. The study con- 
cludes that the galaxies do not trace the underlying dark- 
matter distribution with 61 = 0.93'i'Q Qg and C2 = 1)2/ Bi — 
— 0.3410 ^3 . Observations of the HI bispectrum would allow 
us to carry out such studies at redshifts z > 1, providing fur- 
ther tests of the gravitational instability picture. Further, it 
would be possible to determine the bias of the HI at high 
redshifts, this being a potential probe of galaxy formation. 

The cosmological evolution of the neutral hydrogen 
during the Epoch of Reionization (z = 6 — 20) and the pre- 
Reionization era {z > 20), and the expected redshifted 21 cm 
signal are topics of very intense research (egJSco^^n^^^es 
199(1 iMadau. Meiksin and Reed [l997l: iGnedin fc Ostriker 



1997; Shaver et al."l999'; "Tozzi et al."200d'; 'lliev et al.'"200; 
Furlanetto. Sokasian. fc Hcrng uist 2004; Miral da-Escudi 

2003,1; Zaldarriaga. Furlanetto and Hernquisll 

Loeb fc Zaldarriaea. j2004i ; [BharadwaLfcA^ 



Zaroubi fc Silkl t2004l: IBharadwai and Ali I 
Bharadwai and PandevI l2005t iBarkana fc Loebl l2005£ 
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Barkana fc Loebl l2005bl: iBarkana fc Loebl l2005d: 
Bowmai^Morares fc HewittI I2OO5F ! We note that most 
of these papers have no direct bearing on the cur- 
rent work which is restricted to the Post-Reionization 
Era (2 < 6) where it is reasonable to assume that 
the HI traces the dark matter, and the process of 
gravitational instability has progressed sufficiently to 
produce measurable departure s from the Gaus si an in i- 
tial conditions . iSubramanian fc Pad manabhtirl l)l993fl . 
Kunw^Padmmi^bhai^^^uteam il995fl and 

Bagk^Natl^^PaAnane^heml 1^9^ have considered 
the possibility of detecting HI at z < 6 focusing on the 
possibility of d etecting individual featur es corresponding to 
big HI clumps. iBada and White! ^2003^ have used N-body 
simulations to estimate the two-point statistics of the low z 
HI signal. 

We next present the structure of the paper. Section 2 
reviews the formulae needed to calculate the visibility cor- 
relation signal expected from HI at high z. In Section 3. we 
present the results and discuss their implications. 

Finally we note that, unless stated otherwise, we use 
the values {D„no,^XoMbh'^ , h) = (0.3,0.7,0.02,0.7) for the 
cosmological parameters. Further, we use Q,g = 1 x 10~^ as 
a fiducial value and present results only for this throughout 
the paper. 



2 CALCULATING THE HI VISIBILITY 
CORRELATION SIGNAL 

The quantity measured in radio-interferometric observations 
is the complex visibility V{\J, v) which records only the an- 
gular fluctuations of the speciflc intensity 7^ on the sky. 
This is measured for every pair of antennas in a radio- 
interferometric array. For any pair of antennas at a sepa- 
ration d, we refer to the two dimensional vector U = d/A 
as a baseline. In this paper it has been assumed that all 
the antennas in the array are coplanar, and that the anten- 
nas all point vertically up wards in the direction m with 
m • U = 0. Further, the individual antennas are assumed to 
have a Gaussian beam pattern A{9) — e~* with 
(in radians) i. e. the beam width of the antennas is small, 
and the part of the sky which contributes to the signal can 
be well approximated by a plane. 

Fluctuations in the high redshift HI density con- 
tributes to the visibilities measured at all frequencies be- 
low 1420 MHz. Following BP05 we consider the contribution 
from the HI signal to the correlations 

S2(Ui, Ua, Az.) = {-[/(Ui, V + Ai.)l^(U2, v)) (1) 
and 

S3(Ul,U2,U3,A//i,A!/2) = {V(Ui,I. + Ai/i)V^(U2, 

V + /^1^2)V{\J3,V)) (2) 

expected between the visibilities at different baselines and 
frequencies. Here the two visibility correlation function 
S2(Ui, U2, Ai') denotes the correlation expected between 
the visibilities measured at two different baselines Ui and 
U2 one at the frequency u and the other at v-\-/S.i'. It should 
be noted that although we have shown S2 as an explicit func- 
tions of the frequency differences only, it also depends 



on the central value v which is not shown as an explicit ar- 
gument. The value of v will be clear from the context of the 
discussion. Further, throughout our analysis we assume that 
all frequency differences are much smaller than the central 
frequency ie. Av/u <C 1. The deflnition of the three visibility 
correlation S3 closely follows that of S2. 

It has been shown that the two-visibility cor- 
relation probes the power spect r um of HI fluctu- 
ation s iBharadwai and Sethi 1200 ll . iBharadwai and Ali I 
I 2OO5I) and the three visibility correlation the bispectrum 
iBharadwai and Pandevl2005f) . Considering S2 (Ui , U2 , Ai/) 
first, we note that it has a significant value only when 
I U2 + Ui |< {■k6o)~^ ■ Therefore we restrict the analysis 
to Ui = — U2 = U which we denote as S2(U, Ai/). This is 
related to the redshift space HI power spectrum PHi(k) as 

S2((7,A//) = ^ / dfc||fki(k)cos(fc||r,A//). (3) 
'^'^•^ Jo 

where is the comoving distance to the HI whose 1420 MHz 
radiation is observed at a frequency z^, — dr^/dv, k = 

fclim + (2Vr.)U and /. = 2.5 x 10^ ^ (H^) (M) 

A point to note is that PHi(k) includes the effects of redshift 
space distortions, and is therefore anisotropic along the the 
line of sight m. We further note that 82 is real and isotropic 
le it does not depend on the direction of the baseline U 
whereby we can write S^iU, Au). 

Considering the three visibility correlation 
S3(Ui, U2, U3, Ai/i, Ai/2) next, this has a significant 
value only ff | Ui -I- U2 -I- U3 |< (7r6lo)"^ We therefore 
restrict our analysis to only those combinations of baselines 
for which Ui + U2 + U3 = 0. Further S3 is real and it 
depends only on the shape and size of the triangle which is 
completely specified by the magnitude of the three baselines 
{Ui,U2,U3). We have 

S3(C/l,[/2,?73, Az^l, Az/2) = °4 / rffcl||rffe2|| COs[(fci|| At/1 

LZ TV Tu J 

+ fc2||A!/2)r-l,] X BHi(ki,k2,k3)(4) 

where ki = fci||m-|-(27r/r^)Ui, k2 — fc2||m-(-(27r/ry)U2, 
ks = — (feiii -(-fc2||)m-|- (2-7r/r^)U3 and _BHi(ki, k2, ks) is the 
redshift space HI bispectrum. 

It is next necessary to specify the form of the redshift 
space HI power spectrum Phi and bispectrum Phi- Fluctu- 
ations in the redshifted HI radiation arise from fluctuations 
in the HI density and from HI peculiar velocities, and it is 
necessary to include both these effects in Phi and Phi- We 
assume that the HI traces the dark matter with a possible 
bias, and retain terms upto the quadratic order in the rela- 
tion between the fluctuations in the HI density Sm and the 
fluctuations in the dark matter density 5 

5m = bi5+^S^ (5) 

where bi and 62 are the linear and quadratic bias param- 
eters respectively. A non-zero 62 would indicate nonlinear 
biasing of HI distribution with respect to underlying mass 
distribution. It is also assumed that the peculiar velocities 
are determined by the dark matter which dominates the dy- 
namics. 
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The HI power spectrum PHi(k) is related to the real 
space dark matter power spectrum P{k) as 



PHi(k) = 6?x|i(l + /3^2)2p(fc) 



(6) 



^ being the cosine of the angle between k and the line of 
sight m, and (3 = f{0.m)/bi ~ /bi is the linear dis- 
tortion parameter. The term (1 + in equation Q in- 
corporates the effect of the peculiar velocities, and xui is 
the mean hydrogen neutral fraction. Converting Qgaa to the 
mean neutral fraction xm = pni/pH = Qgas/^b gives us 
xm = 50ngash'^{0.02/ni,h'^) or = 2.45 x 10"^. It should 
be noted that | 5 |^ 1 and terms higher than have been 
dropped in the power spectrum. 

For the HI bispectrum we assume 



PHi(ki,k2,k3 ,z) = a;Hii3s(ki,k2,k3 ,z) 



(7) 



whe re is the redshift space bispectrum calculated by 
199 dVerde eraP 1998 1998. We have used the fo rm of Bs 
given m equation (11) of ( 199 dVerde et al. 19'93) . the ex- 
pression being quite lengthy we do not reproduce it here. 
We note that Ps(ki,k2,k3) has a quadratic dependence on 
the dark matter power spectrum P(fc) , it depends on the bias 
parameters hi and 62, and on the parameters of the back- 
ground cosmological model. The value of B3 also depends 
on the shape of the triangle formed by (ki, k2, ks). 

The fact that the neutral hydrogen is in discrete clouds 
makes a contribution which we do not include here. This 
effect originates from the fact that the HI emission line 
from individual clouds has a finite width, and the visibility 
correlation is enhanced when Ai/ is smaller than the line- 
width of the emission from the individual clouds. Another 
important effect not included here is that the fluctuations 
become significantly non-linear at low z. Both these effects 
have be en studied for the power spectrum using simulations 
tBharadwai fc Srika nt 2004") . 

The simple analytic treatment adopted here suffices for 
the purposes of this paper where the main focus is to esti- 
mate the magnitude and the nature of the HI signal, and 
to investigate the feasibility of using such observations to 
probe structure formation at high redshifts. 



3 RESULTS AND DISCUSSION 

We present results for the two and three visibility correla- 
tion signal expected from HI at high z. In calculating the 
expected signal we have used the GMRT parameters, but 
it is straight forward to scale the results presented here 
to make the visibility correlation predictions for other ra- 
dio telescopes. The only telescope parameter which enters 
equations J^l and Q is 60 , which is the beam size of the 
individual antenna in the array. Further, it should be noted 
that 60 ~ 0.6 X 6fwhm . The value of 60 will depend on the 
physical dimensions of the antennas and the wavelength of 
observation. For the GMRT 6I0 = 1° at 325 MHz. We scale 
this using So oc A to obtain 60 corresponding to different 
observation frequencies. We note that the smallest baseline 
at the GMRT is a little smaUer than U = 100 while the 
largest is around U ~ 10, 000, and we present results for U 
in the range 10 — 10, 000. It may be noted that unless men- 
tioned otherwise ,we use the values (61 , 62 )~ (1.0 , 0.5) for 
the bias parameters. 




Figure 1. This shows the expected two visibility correlation sig- 
nal [S2(f/, 0)]^/'^ as a function of U for the different values of the 
central frequency shown in the figure. 



Figure Q shows the [S2{U,0)]^^^ , the two visibility cor- 
relation signal expected when the two visibilities are at the 
same baseline and frequency. The signal is strongest at small 
baselines, and it is around 0.1 mjy at 325 MHz ie. {z = 3.4). 
The signal falls with increasing U, and the U dependence 
refiects the shape of the power spectrum. Further, the sig- 
nal increases with increasing central frequency or decreasing 
redshift,. This is a consequence of the fact that the matter 
power spectrum grows with time. 

To estimate the magnitude of the expected three visi- 
bility correlation signal S3((7i, U2, U3, Aui, Ai^2), we first re- 
strict our analysis to equilateral triangles for which the size 
of the baseline U completely specifies the triangle, We fur- 
ther restrict our analysis to Ai^i = Ai/2 ~ Av in which case 

[S3(C/,0)]^''^ 
[S2(C/,0)ll/2 



we can use the notation S3(f/, Ai^). The ratio 
(Figurel^J quantifies the relative strength of S3 as a fraction 
of S2. The first point to note is that the three visibility corre- 
lation signal is comparable in magnitude to the two visibility 
correlation signal. The ratio of the signal strengths is of order 
unity for the small baselines {U < 1000) at 325 MHz and it 

is somewhat higher at higher frequencies. Further, S3'''' falls 
1/2 

faster than S2 with increasing U. It is possible to under- 
stand both these features of S3 in terms of the fact that the 
bispectrum is quadratic in the power spectrum whereby it 
grows faster with time and it has a steeper k dependence at 
large k. 

The bispectrum Ps(ki,k2,k3) depends on the shape 
and size of the triangle formed by k2 , k2 and ks . This de- 
pendence varies with the bias parameters (Bi,62), and this 
can be used to observationally determine the values of the 
bias parameters. To study how the shape and bias depen- 
dence manifests itself in S3, we have considered bilateral 
triangles for which two sides are of length U with opening 
angle 120° and the third side of length \^U. We study the 
ratio [SsiU, Aiy)]l^^ /[SaiU, Aiy)]l^^ for Au = where the 
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Figure 2. This shows the ratio jg^|y'Q|ji/2 ^ function of U for 
the different values of the central frequency shown in the figure. 
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Figure 3. This shows the ratio [SsiU, Au)]]^^ /[SaiU, Au)]]^^ for 
Au = 0, where E refers to equilateral triangles and B to bilateral 
triangles. Two sides of the bilateral triangle are U (at 120°) and 
the third is \/3C^- The central frequency is 610MHz. The results 
are for the different combinations of the bias parameters (f)i,f)2) 
shown in the figure. 



Figure 4. This shows the expected two visibility correlation sig- 
nal [S2{U, Au)]^/'^ as a function of Au for the different values of 
U shown in the figure at a central frequency of 610 MHz 
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subscripts E and B refer to equilateral and bilateral triangles 
respectively. Figure|21shows the results at 610MHz {z — 1.3), 
the behaviour at other frequencies is similar. We find that 
S3 does not show any significant dependence on the triangle 
shape at baselines smaller than U ~ 1000. There is a shape 
dependence at larger baselines where S3 decreases faster for 
the equilateral triangles than for the bilateral triangles. Fur- 
ther, the ratio has a distinct dependence on the values of the 
bias parameters. The first point to note is that the ratio is 
independent of 61 for 62 ~ 0. For other values of 62 we find 
that the U dependence of the ratio [S3 ((7, 0%^^ /[SaiU, 0)]]^^ 
at the large baselines is sensitive to both bi and 62. It should, 
in principle, be possible to use observations of the shape de- 
pendence of S3 to determine the bias parameters. 

We next consider the correlation between visibilities at 
different frequencies. The results are presented at fixed val- 
ues of the baselines, and Figures|l]and|^show [S2(U, Av)]^^'^ 
and [53(11, Ai/)Y^^ respectively, both at the central fre- 
quency of 610 MHz. We find that the correlations decay 
rapidly as the frequency separation Av is increased. We have 
S{U,Siy) = at Ai/ ~ 1.3MHz for U = 100, and there is a 
weak anticorrelation S2 < for larger Av. The visibilities 
at larger baselines get decorrelated faster with increasing U, 
and the zero crossing is at ~ 0.7 MHz at ?7 = 100. The three 
visibility correlation shows a similar behavior, the main dif- 
ference being that [S3(I7, Ai^)]^'''' decays much faster than 



Figure 5. This shows the expected three visibility correlation 
SIB nal [S3(C/, Ai/)]i/3 

as a function of Au for the different values 
of U shown in the figure at a central frequency of 610 MHz 



[S2{U, Au)]^^^ with increasing Au. It may be noted that the 
overall behavior is similar at other values of the central fre- 
quency, with the difference that S2 and S3 decay faster with 
increasing Av at lower values of the central frequency. 

An important issue which we have not addressed here 
in any detail is the extraction of the HI signal from various 
contaminants which are expected to dominate at low fre- 
quencies. The issue of extracting the red shifted 21 cm signal 
has received considera b le attentio n (eg. |ShaveiLet_alj[l999|; 
i DiMatteo. et alJl2004l : lOh et alJ l2003llMorales and Hewitt 



(2004 

2004 



IZaldarria.a , Furl anetto and HernquistI 

ICoorav fc Furlanetto 20(53 : ISantos et al 
IWang ■ Teeinark .Santos fc Knoxl l2005l 



20051 ; 

Morales .Bowman fc Hewitt! l200,4 l 



The contribution 
from the continuum sources of contamination at two 
different frequencies is expected to be correlated even at 
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large l\v, whereas the HI signal is found to be uncorrelated 
beyond something like 1.3 MHz or even less depending on 
the value of U . It is then, in principle, straightforward to 
fit the visibility correlations S2 and S3 at large /S.v and 
remove any slowly varying component thereby separating 
the contaminants from the HI signal. We also use this 
opportunity to note that this is a major advantage of 
using visibility correlations as compared to the angular 
power spectrum which has been found to have substantial 
correlations even at two frequencies separated by ~ 10 MHz 
JSantos et alj|2005h . 

System noise is inherent in all radio interferometric ob- 
servations. System noise and the observation time needed 
to detect the HI signal is another important issue. Here we 
briefly discuss the signal to noise ratio and approximate in- 
tegration time required to detect the HI signal for S2 and 
S3. We consider an array of A'^ antennas, the observations 
lasting a time duration t, with frequency channels of width 
5v spanning a total bandwidth B. It should be noted that 
the effect of a finite channel width 5p has not been included 
in our calculation which assumes infinite frequency resolu- 
tion. This effect can be easily included by convolving our 
results for the visibility correlation with the frequency re- 
sponse function of a single channel. Preferably, 5v should 
be much smaller than the frequency separation at which 
the visibility correlation become uncorrelated. We use <S 
denote the frequency separation within which the visibili- 
ties are correlated, and beyond which they become uncor- 
related. We use N2 and N3 to denote the rms. noise in S2 
and S3, we have assumed that thermal noise are Gaussian 

random fields. It is well known that N2 = 377 

ijThompson. Moran fc SwensonlFlflSfif l. and we have N3 ~ 

where Tsys is the system tempera- 



{il^t)3/2 ' 



ture and A^j is the effective area of a single antenna.The 
noise contribution will be reduced by a factor 1 / No if we 
combine No independent samples of the visibility correla- 
tion. A possible observational strategy for a preliminary de- 
tection of the HI signal would be to combine the visibility 
correlations at all baselines and frequency separations where 
there is a reasonable amount of signal. This gives A'^o = 
[N{N - l)/2] {B/5v){S/Sv) for the two visibility correla- 
tion and No = [N{N - 1){N - 2)/6] {B/Su) {S/5uf for the 
three visibility correlations. Combining all of this we have 

''^'^ and [N3]^/3 _ 
Using values for the 



[j^^jl/2 _ (2i^^) 

\ ^'f ) 



N(N-1)BS 
1/6 



N(N-1)(N-'2)B 52 



1 

172 • 



GMRT at ~ 610MHz, {2kBTsYS /Aef) = 144 Jy, B = 
16 MHz, there are 14 antennas within U < 1000 where the 
signal is strong and 5 = 0.5 MHz beyond which the signal is 
uncorrelated we find that it is possible to achieve noise lev- 
els of [N2]^''^ ~ 0.03 mJy and [N3]^/^ ~ 0.045 mjy which 
are below the signal with 200 hrs of integration. Our present 
estimations indicate that 200 to 1000 hours of observation 
are needed to detect the HI signal. 
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